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Abstract. We provide a general macrostatistical formulation of nonequilibrium steady 
states of reservoir driven quantum systems. This formulation is centred on the large scale 
properties of the locally conserved hydro dynamical observables, and our basic physical 
assumptions comprise (a) a chaoticity hypothesis for the nonconserved currents carried 
by these observables, (b) an extension of Onsager's regression hypothesis to fluctuations 
about nonequilibrium states, and (c) a certain mesoscopic local equilibrium hypothesis. 
On this basis we obtain a picture wherein the fluctuations of the hydro dynamical variables 
about a nonequilibrium steady state execute a Gaussian Markov process of a generalized 
Onsager-Machlup type, which is completely determined by the position dependent trans- 
port coefficients and the equilibrium entropy function of the system. This picture reveals 
that the transport coefficients satisfy a generalized form of the Onsager reciprocity relations 
in the nonequilibrium situation and that the spatial correlations of the hydro dynamical 
observables are generically of long range. This last result constitutes a model-independent 
quantum mechanical generalization of that obtained for special classical stochastic systems 
and marks a striking difference between the steady nonequilibrium and equilibrium states, 
since it is only at critical points that the latter carry long range correlations. 
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1. Introduction 

The statistical thermodynamics of nonequilibrium steady states or, more generally, dy- 
namically stable ones, of reservoir driven macroscopic systems 6 is a key area of the natural 
sciences, with ramifications for condensed matter physics [1-4], chemistry [5] and biology 
[6]. At the phenomenological and heuristic levels, there is an abundant literature on this 
subject. At the level of mathematical physics, however, the subject is still at an exploratory 
stage. In the classical regime, two types of rigorous approaches have been made to it. The 
first is centred on the hypotheses that the macroscopic properties of complex systems are 
yielded by the model of classical Anosov dynamical systems [7,8]. This hypothesis is de- 
signed to capture the chaoticity that underlies macroscopic irreversibility, and it has been 
shown to lead to nonequilibrium generalizations both of the Onsager reciprocity relations 
[8] and of the fluctuation-dissipation theorem [7] . A second approach is centred on micro- 
scopic treatments of stochastic (non-Hamiltonian) dynamical models [9-11], which are also 
designed to capture the chaoticity underlying macroscopic irreversibility. The treatment 
of these models has led to some interesting developments, and Ref. [11] has provided a 
dynamically based picture of the hydro dynamical fluctuations about their nonequilibrium 
steady states. Moreover, in the case of a certain particular model, namely the symmetric 
exclusion process, it has been shown that the nonequilibrium steady state has long range 
density correlations [9-11] and that the probability distribution of its large scale density 
field is determined by an explicitly specified and highly nontrivial nonequilibrium gener- 
alization of its free energy [10, 11]. In the quantum regime, a natural dynamically based 
definition of nonequilibrium steady states of reservoir driven systems has been formulated 
[12, 13] at the microscopic level. 

In the present article we set out a different approach to the subject, which is quan- 
tum macro statistical in that is is centred on the hydro dynamical observables of reservoir 
driven quantum systems. This approach, which was briefly sketched in Ref. [14], parallels 
the one we have previously made to the nonequilibrium thermodynamics of conservative 
quantum systems [15, 16], where it yielded an extension of the Onsager reciprocity re- 
lations to a nonlinear regime. In general, the quantum macrostatistics is designed, like 
Onsager's [17] irreversible thermodynamics and Landau's fluctuating hydrodynamics [18], 
to form a bridge between the microscopic and macroscopic pictures of matter, rather than 
a deduction of the latter from the former. Indeed, accepting Boltzmann's hypothesis of 
molecular chaos [19], we take the view that such a derivation is not even feasible for real- 
istically interacting systems, since this chaos renders the microscopic equations of motion 
intractable over periods substantially longer than the intervals between successive colli- 
sions . Thus, the microscopic equations of motion must necessarily be supplemented by 
further assumptions in order to interconnect the quantum and phenomenological properties 

b A very simple example of such a state is the stationary one of a solid rod, whose ends 
are coupled to thermostats of different temperatures. 

c This view is supported by the fact that the rigorous derivations of Boltzmann equa- 
tions from the Hamiltonian dynamics of both classical [20] and quantum [21] systems are 
applicable only over microscopic times of the order of the interval between successive col- 
lisions of a particle. For longer times, the chaos bars the way to further analysis of the 
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of matter. In fact, the key physical assumptions of our macrostatistical project concern 
only very general, model-independent properties of many-particle systems. Specifically, 
they comprise 

(A) an extension of Onsager's regression hypothesis [17], to the effect that the hydrodynam- 
ical fluctuations about nonequilibrium steady states are governed by the same dynamical 
laws as the 'small' perturbations of the hydrodynamical variables about their steady values; 

(B) a certain mesoscopic local equilibrium hypothesis; and 

(C) a chaoticity hypothesis for the nonconserved currents carried by the locally conserved 
hydrodynamical observables. 

These assumptions may be regarded as the 'axioms' of our theory. The physical consider- 
ations that underlie them will be discussed, along with their formulation, in Sections 4.1, 
4.1 and 4.4. In fact, the hypothesis (C), like Boltzmann's Stosszahlansatz and its subse- 
quent developments [7-11], exploits the consequence of the very chaos that obstructs the 
analytical dynamics of realistically interacting many-particle systems. 

The principal results that we obtain by supplementing the Schroedinger dynamics of 
many-particle systems by the 'axioms' (A)-(C), together with certain technical assump- 
tions, are the following ones (I)- (III), which we claim to be new, at least on the level of a 
rigorous, general, model-independent quantum theory of nonequilibrium steady states. 

(I) The spatial correlations of the hydrodynamical observables are generically of long range. 
This comprises a quantum mechanical generalization of that obtained from both rigorous 
microscopic treatments of certain classical stochastic models [9-11] and from heuristic 
treatments [23, 24] of Landau's fluctuating hydrodynamics. Most importantly, it marks 
a qualitative difference between equilibrium and nonequilibrium steady states, since the 
hydrodynamical correlations in the former states are generically of short range, except at 
critical points. 

(II) The transport coefficients satisfy a generalized, position-dependent version of the On- 
sager reciprocity relations. Thus, this result extends Onsager's irreversible thermodynam- 
ics from the neighbourhood of equilibrium to that of nonequilibrium steady states. 

(III) The hydrodynamical fluctuations execute a classical Gaussian Markov process, of a 
generalized Onsager-Machlup (OM) type [22]. Thus this result extends the OM theory 
from the regime of fluctuations about thermal equilibrium to that of fluctuations about 
nonequilibrium steady states. A similar result was obtained for certain classical stochastic 
models in Ref. [11]. 

Let us now briefly describe the macrostatistical strategy we employ to obtain these 
results. We take our model to be an A^-particle quantum system, E, that is confined to a 
bounded open connected region, fi^r, of a <i-dimensional Euclidean space, X, and coupled 
at its boundary, dVt at, to an array, TZ, of quantum mechanical reservoirs. E is thus an open 
system, while the composite (E + 1Z) is a conservative one. Since we shall have occasion 



microscopic equations of motion. 
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to pass to thermodynamic and hydrodynamic limits where its particle number tends to 
infinity, we take N to be a variable parameter of the system. We assume that its particle 
number density v := N/Vo\(Qn) is ^-independent and that Qn is the dilation by a factor 
Ln of a fixed, N- independent region of unit volume. Thus Qn = L^Q := {Lnx\x&Q} 
and 

Ln = (iV» Vd - (1-1) 

For the hydrodynamic description of E, we take Ln to be the unit of length. Thus, Q is 
the region occupied by the system in the hydrodynamical picture. 

We assume that, in that picture, £ evolves according to a phenomenological law gov- 
erning the evolution of a set of locally conserved classical fields qt(x) = (qi jt (x), . ., q m ,t(x)), 
which correspond to the densities at position x and time t of the extensive thermody- 
namic variables'^ of the system. We denote the associated currents of qt(x) by jt(x) = 
(ji,t(x), ■ -,jm,t( x ))- Thus, q t satisfies the local conservation law 

^+V.j,(x) = 0. (1.2) 

We assume that its phenomenological dynamics is governed by a constitutive equation of 
the form 

jt(x) = J(q t ;x), (1.3) 

where J is a functional of the field q t and the position x. Thus, by Eqs. (1.2) and (1.3), 
q t evolves according to an autonomous law 

^^=^(q t ;x):=-V.J(q t ;x), (1.4) 

subject to boundary conditions determined by the reservoirs. We assume that this phe- 
nomenological law is invariant under scale transformations x^\x, t^X k t for some constant 
k. A simple example for which this assumption is valid, with k = 2, is that of nonlinear 
diffusions, where J takes the form 

J(q t] x) = -K(q t (x))Vq t (x), (1.5) 

K being an m-bj-m matrix [-ftTfcz], which acts by standard matrix multiplication on Vqt- 
In this case, the phenomenological equation (1.4) takes the form 

^ = V.(K(q t )Vq t ). (1.6) 

We shall base some of our explicit calculations on this case and, in particular, we shall 
henceforth assume that the scaling exponent k is equal to 2. A simple consequence of this 



d We provide a characterization of these variables in Section 2.2 along lines previously 
formulated in Ref. [15]. 
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assumption is that, since Ln is the unit of length for the hydrodynamical picture, L 2 N is 
the unit of time for this picture. 

We assume that, in general, the dynamics described by Eq. (1.4) is dissipative, in that 
the m-component field qt(x) relaxes eventually to a unique time-independent form q(x), 
which thus corresponds to a steady hydrodynamical state. By Eq. (1.3), the corresponding 
steady m-component current, j(x), is then J~(q;x). 

By Eq. (1.4), the linearised equation of motion for 'small' perturbations, 8q t (x), of 
q(x) is simply 

— Sq t (x) = CSq t (x) := q^^(q + *8q t ;x)\ x=0 , (1.7) 
while, by Eq. (1.3), the corresponding increment in the m- component current j(x) is 

d 

Sj t {x) =KSq t {x) := -^J(q + XSq t ;x) lx =o- (1-8) 

We note that, by Eqs. (1.4), (1.7) and (1.8), 

C=-VX. (1.9) 

Further, in the case of nonlinear diffusions, it follows from the identification of the r.h.s.'s 
of Eqs. (1.4) and (1.6) that Eq. (1.7) yields the following formal equation for C. 

[C X ](x) = V.(ir( 9 (x))Vx(x) + [K'(q(x)) X (x)} Vg(x)), (1.10) 

where x is a single column matrix function of position and K'(q) is the derivative of K(q), 
i.e. its gradient with respect to q: thus [K'(q)x(q)]ki = YlT=i[ d ^kl(q)/dq r ]Xr(q)- 

In order to relate the phenomenological dynamics given by Eqs. (1.4) and (1.7) to the 
underlying microscopic quantum mechanics of E, we assume that qt(x) is the expectation 
value of a set of locally conserved quantum fields qt(x) = (qi it (x), . ., q m ,t(x)) as rescaled 
for the hydrodynamical picture and in a limit in which N, and hence Lat, becomes infinite. 
Correspondingly, we formulate the fluctuations £t(x) of this m-component quantum field 
qt(x) about its mean on the same macroscopic scale and with a standard normalization, 
subject to the above-described assumptions (A)-(C). 

On this basis, we establish that & executes a Gaussian Markov process represented 
by a generalized Langevin equation of the form 

J^t(s) = + &*(*). (l-H) 

where bt(x) is a white noise whose autocorrelation function is of zero range with respect 
to position as well as time. Thus, £t executes a generalized Onsager-Machlup process. 
We employ this result to infer that the spatial correlations of the fluctuation field £ in 
nonequilibrium steady states are generically of long range. In this way we derive the above 
results (I)- (III) from our basic macrostatistical assumptions. 
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We present our treatment as follows. In Section 2 we formulate the quantum statistical 
thermodynamical model of the composite system (H + 1Z) at both microscopic and macro- 
scopic levels. This formulation provides general specifications of the nonequilibrium steady 
states of the model and also of the locally conserved quantum fields qt and associated cur- 
rents jt pertinent to its hydrodynamic description. Here, in accordance with the general 
requirements of quantum field theory [25], we assume that these are distribution- valued op- 
erators. In Section 3 we relate the classical hydrodynamical variables, q t and j t , and their 
fluctuations, £ t and r\ t , about a nonequilibrium steady state to these quantum fields and 
currents; and we obtain sufficient conditions for the fluctuations £t to execute a classical 
stochastic process. In Section 4 we formulate our regression and local equilibrium hy- 
potheses for this process and note that these, together with the assumption of microscopic 
reversibility for the composite (E+7£), yields a canonical extension of Onsager's reciprocity 
relations to the nonlinear hydrodynamical regime. In Section 5 we extend our local equi- 
librium hypothesis to the fluctuating currents, rjt, and formulate our chaoticity hypothesis 
for these currents. We then establish that the assumptions of the regression hypothesis, 
local equilibrium and chaoticity imply the field £ t executes a generalized Onsager-Machlup 
process represented formally by Eq. (1.11). In Section 6 we obtain an explicit formula for 
the two-point function for this process in terms of the equilibrium entropy density function 
and the transport coefficients of the system, and we infer therefrom that the static corre- 
lations of the hydrodynamical fluctuation field £ are generically of non-zero range on the 
macroscopic scale and hence of long (infinite!) range on the microscopic one. We conclude 
in Section 7 with some general observations about the results of this article and of their 
possible generalizations to less restrictive conditions than those assumed here. We leave 
the proofs of some technical Propositions to four Appendices. 

2. The Quantum Model. 

We take our model to be the open quantum system, E, briefly described in Section 1. 
Thus, E is a system of N particles, which occupies a bounded open connected region, Ojv, 
of a (i-dimensional Euclidean space X and is coupled at its surface, <9Qat, to an array, 1Z, 
of reservoirs. Here Qn is the dilation by a factor Ln of a region, O, of unit volume and 
Ln is given by Eq. (1), which represents the iV-independence of the particle density of E. 
We assume that the composite quantum system E( c ) := (E + 1Z) is conservative and that 
all its interactions are invariant under spatial translations and rotations. 

2.1. The Microscopic Picture. We formulate this picture in standard operator 
algebraic terms, denoting the C*-algebras of bounded observables of E and E( c ) by A and 
B , respectively. We assume that A is a subalgebra of B and that it is isomorphic to the 
VF*-algebra of bounded operators in a separable Hilbert space 7i, which comprises the 
square integrable functions f(xi, . ,xn',si,. ., s^) (appropriately symmetrized or antisym- 
metrized) of the positions {xj} and the spins {sj (= ±1)} of its particles. The unbounded 
observables of E are represented by the unbounded self-adjoint operators affiliated to A, 
i.e. by those whose spectral projectors belong to this algebra. The states of this system are 
represented by the density matrices in Ti, and the expectation value of an observable, A, 
of E for the state p is Tr(pA). In general we denote this expectation value by p(A)=(p; A), 
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and we employ the corresponding notation for £( c ). 

The Wigner time reversal operator, which serves to reverse the velocities and spins of 
the particles of E, is defined to be the antilinear transformation of 7i given by the formula 

(Tf)(x 1 , . , x N ; si, . ., s N ) = f(x 1 , . , x N ; -si, . ., -s N ) V/eW, (2.1) 

where the bar denotes complex conjugation. Thus, T implements an antiautomorphism 
T4 of A, defined by the formula 

t a A = TA*T V AeA. (2.2) 

We assume that the dynamics of the composite system £( c ) is given by a one-parameter 
group, {o^|tGR} := a(R), of automorphisms of B. Further, we assume that this dynamics 
is reversible, i.e. that B is equipped with an antiautomorphism r, which reduces to T4 on 
A and implements time reversals according to the prescription 

ra t T = a-f (2.3) 

The evolution of the observables of £ is given by the isomorphisms of A into B obtained 
by the restriction of a(R) to the former algebra. 

2.2. Thermodynamic Variables and Potentials. In order to formulate the 
thermodynamic observables and potentials of £ we pass, for the moment, to the situation 
where it is decoupled from the reservoirs 1Z and thus becomes a conservative system, whose 
dynamics is given by a one-parameter group, {a[ ^ |tGR}, of automorphisms of A. In this 
situation, its canonical equilibrium state, p, at inverse temperature (3 is characterized by 
the Kubo-Martin-Schwinger (KMS) condition [26] 

(p; [a ( t 0) A l ]A 2 ) = (p; A 2 a ( t % h0 A 1 ) V A u A 2 e A; t e R. (2.4) 

Most importantly, this condition survives the thermodynamic limit where iV tends to 
infinity and the particle density v remains finite [26]. Moreover, in this limit 6 , the system 
may support different states that satisfy the condition. The set of these states is convex, 
and its extremal elements may naturally be interpreted as the pure equilibrium phases for 
the inverse temperature f3 [15, 29]. 

We assume that £ has a linearly independent set of extensive conserved observables 
Q = (Qi, • -,<5n), which intercommute-f up to surface effects and satisfy the following 

e The model of the infinite system is formulated, in a standard way, in terms of its 
C*-algebra of quasi-local bounded observables [15, 26-28]. Its states are then positive 
normalized linear functionals on that algebra. 

f The assumption of intercommutativity is not universally fulfilled. It is violated, for 
example in the case where Qk and Qi, say, are different components of the magnetic 
moment of £. In such cases, some aspects of our treatment would have to be refined. 
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condition of thermodynamical completeness [15]:- in the limit N-^oo, the pure phases are 
labelled by, i.e are in one-to-one correspondence with, the expectation values q\,. .,q m of 
the global densities of Qi, ■ -,Qm, respectively. The resultant set of classical, intensive 
thermodynamical variables of E is then q = (qi, . ., q m ). In general, we take Q± to be the 
Hamiltonian of the system: correspondingly, q\ is its energy density. 

The equilibrium entropy density, in the limit iV^oo, is a function, s, of q, which 
may be formulated by standard methods of quantum statistical mechanics [15, 27]. The 
classical equilibrium thermodynamics of the system is then governed by the form of s(q). 
The demand of thermodynamical stability ensures that this function is concave. We define 
the thermodynamic conjugate of q^ to be 9k = ds(q)/dqk- Thus, denoting the element 
(#!,.., # m ) of R™ by 

9 = s'(q), (2.5) 

the derivative of s(q), i.e. its gradient in (/-space. Correspondingly, the second derivative, 
s"(q), of this function is the Hessian [d 2 s(q)/dqkdqi]. We assume throughout this treatment 
that the system is in a single phase region, i.e. one where s is infinitely differentiable, 
where the function q^9(q) is invertible and where, for each value of q, the matrix s"(q) is 
invertible. We define 

J(q) := -s'\q)-\ (2.6) 
which, in view of the concavity of s, is a positive matrix. 

2.3. The Reservoir System 1Z. We assume that 1Z comprises a set, {TZj}, of 
spatially disjoint reservoirs, such that each IZj is placed in contact with a subregion dQ^^j 
of <9fi n and [J jOQn^j = dQpj. Further, we assume that each TZj has a ther mo dynamically 
complete set of global extensive conserved observables (Qj,i, ■ ., Qj, m ) that are the natural 
counterparts of Qi, . ., Q m , respectively, in that, when E and TZj are placed in contact, the 
observables (Qk + Qj,k) of E^ c ^ are still conserved. Correspondingly, the thermodynamic 
control variables of IZj conjugate to Qj are the same as those of E, namely 9. We denote 
by u>j(9j) the equilibrium state of IZj for which its 9- value is 9j. 

2.4. Nonequilibrium Steady States of E^ c ^ . Returning now to the situation where 
E is an open system, we assume that this is prepared according to the following prescrip- 
tion. E and the reservoirs {TZj} are independently prepared in the remote past in states po 
and {u)j(9j)}, respectively, where po is normal and the values of 9j generally varies from 
reservoir to reservoir: thus, in general, the reservoirs {TZj} are not in equilibrium with one 
another. Following this preparation the systems E and TZ are then coupled together and 
the resultant conservative composite evolves freely according to the dynamics governed by 
the automorphisms ck(R). We assume that, as established under suitable asymptotically 
abelian conditions [12, 13], this dynamics acts so as to drive the system 3 E^ c ^ into a ter- 
minal po-i n d e P en <ient state 0(= w* — \\va t ^ 00 a^[pQ®jujj{9j)} ), whose restriction to A is 



9 The same result has been also obtained constructively [30] for certain models, which 
however are too rudimentary for our present purposes. In particular, the version of E there 
is just an multi- level atom. 
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normal. This state is uniquely determined by the states {uj(9j)}. Accordingly, we take 
4> to be the nonequilibrium steady state of T,^ stemming from the specified preparation, 
and we denote its GNS triple by (7^, 7r, <E>). 

We note that, in view of the stationarity of 0, the automorphisms a(R) are imple- 
mented by a unitary representation U of R in according to the prescription [31] 

n(a t B) = U t n{B)U^ V BeB, teR, (2.7) 

where U is defined by the formula 

U t n(B)$ = n(a t B)$ V BeB, teR. (2.8) 

Since Eq. (2.8) is applicable to the subalgebra A of £>, the dynamics of the open system 
E, in the normal folium of 0, is given by the isomorphisms implemented by U of tt(A) into 
tt(B). 

Moreover, this prescription extends to the unbounded observables of E for the fol- 
lowing reasons. Since the restriction of <p to A is normal, so too, by Eq. (2.7), are the 
representations 7r and 7roct t . It follows [32] that these representations have canonical ex- 
tensions to the unbounded observables, S, of E according to the prescription that, if {E\} 
is the family of spectral projectors of S, then those of n(S) and n(a t S) are {rc(E\)} and 
{n(a t E\) = U t ir(E\)Uf }, respectively. Hence, the extension of the formula (2.7) to the 
unbounded observables takes the form 

n(a t S) = U t n(S)U- t (2.9) 

for all unbounded observables S of E. 

2.5. The Fields q and the Currents j. We assume that, in the GNS representation 
tv for the nonequilibrium steady state 0, the m-component extensive thermodynamical 
observable Q has a position-dependent, locally conserved density q(x) — (qi(x), . ., g m (x)), 
with associated current density j(x) = (j\(x),. .,j m (x)). Thus the q^s and j^s are 
quantum fields and, in accordance with the general requirements of quantum field theory 
[25], we assume that they are distributions' 1 , in the sense of L. Schwartz [33]. 

We formulate these distributions in terms of the Schwartz spaces, V(Qn) and XV (^jv)j 
of real, infinitely differentiable scalar and R d -vector valued functions, respectively, on X 
with support in Q N . We define V m (Q N ) and T>y(Q N ), respectively, to be the real vector 
spaces given by their m'th topological powers, equipped with the operations of binary 
addition and multiplication by real numbers given by the formula 

A(/l, • fm) + A'(/i, • , f'J = (Afl + X'fi ■ , Xf m + A'/m) 

In concrete cases, it is a simple matter to verify that the explicit formulae for these 
fields and currents are indeed distributions. For example, the number density operator at 
position x is simply ^2^ =1 S(x — x r ), where x r is the position of the r'th particle. 
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V A, A'eR, f k , f' k eV(Q) or V V (Q), k = 1, . , m. 

We denote by £> /m (fi7v) and V^^n) the topological dual vector spaces of £> m (fiAr) and 
Vy (O^) respectively. Evidently, these are spaces of distributions (cf. [33])). 

We assume that the m-component fields q(x) and j(x) are operator valued elements of 
"D /m (07v) and Py^fijv), respectively. For simplicitly, we also assume that the components, 
of g are invariant under time- reversals 1 , i.e that they commute with the Wigner time 
reversal operator T. 

The algebraic properties of the field q(x) are governed by the forms of the commutators 
[qk(x),qi(y)]- We assume that these take the following form, which is readily verified by the 
use of standard formulae in the case where q± is the energy density of the system and the 
remaining q k s are the particle number densities for the different species of its constituent 
particles. 

Em 
c k i r j r (x).VS(x - y), (2.10) 
r—l 

where the c's are iV-independent constants. This formula evidently accords with our 
assumption that Q k s intercommute, up to surface effects: indeed it implies that their 
commutators are just the integrals of currents over OQn- 

We denote by q(f) and j(g) the 'smeared fields' obtained by integrating the distribu- 
tions q and j against test functions / = (/i, . ., f m ) and g = (g\, . ., g m ), which belong to 
the spaces V m (Q N ) and Vy(Q N ) respectively. Thus 

= E fc=1 / dxq k {x)f k {x) (2.11) 

and 

^9) = J2 k=1 / dxj k {x).g k {x). (2.12) 

In general, these smeared fields are unbounded observables, affiliated to the algebra A. 
Therefore, by Eq. (2.7), their evolutes at time t, which we denote by qt(f) and jt(g), are 
their transforms implemented by the unitary operator Ut- Thus, they are the smeared 
fields corresponding to distribution valued operators qt(x) = Utq^U^ 1 and jt(x) = 
U t j{x)U^ 1 , respectively; and the analogous statement may evidently be made for their 
components q k ,t(x) and jk,t( x )- F° r notational convenience, we shall sometimes denote 
Qt(x), q t (f), jt(x) and j t (g) by q(x,t), q(f,t), j(x,t) and j(g,t), respectively. 

We assume that the cyclic vector $ for the state lies in the domain of all monomials in 
the smeared fields qt(f) and jt'(g) and that the resultant vector values of these monomials 
are continuous in the /'s, g's, t's and £"s. 



1 Standard examples of time- reversal invariant q k s are the local number and energy 
densities of many-particle systems. 



11 



Since j is the current associated with q, the local conservation laws for the latter field 
may be expressed in the form 

Qt(f) ~ Uf) = [ duj u {Vf) VMGRJG V m (Q N ). (2.13) 



2.6. The Hydrodynamical Scaling. We assume that the hydrodynamical ob- 
servables of the open system E comprise just the m-component field q, as viewed on the 
scale where the unit of length is L^. Thus, on this scale, the system is confined to the 
fixed region O. Further, in accordance with our assumption, following Eq. (1.6), that the 
macroscopic dynamics is invariant under space-time scale transformations x^Xx, t— >A 2 t, 
we assume that L 2 N is the unit of time corresponding to the length unit L^. Hence, in the 
normal folium of the nonequilibrium steady state cf>, the m-component hydrodynamic field 
is represented by the distribution valued operator 

q t (x) :=q(L N x,L 2 N t). (2.14) 

It follows from this equation and Eq. (2.11) that the smeared hydrodynamic field obtained 
by integrating qt(x) against a D m (fi)-class test function / is 

&(/) = q(f {N) ,L 2 N t), V feV m (Q), teR, (2.15) 

where f( N ' (G£> m (fiAr)) is related to / according to the formula 

/ (x) = L^ d f(L^x) V xen N . (2.16) 

Since the scale transformation (x, t)—>(LNX, L 2 N i) sends q to q, it follows that the local 
conservation law (2.13), or formally dq t (x)/dt = — V.jt(x), will be preserved if its sends 
j t (x) to ] t (x), where 

jt(x) := L N j(L N x,L 2 N t). (2.17) 

It follows from this formula and Eq. (2.12) that the smeared field obtained by integrating 
jt(x) against a Vy(Q)- class test function g is 

Ug) = Kg {N \L 2 N ti (2.i8) 

where 

g( N \x) = L 1 - d g(L- 1 x). (2.19). 

In view of Eqs. (2.15) and (2.18), it is a simple matter to confirm that the local conservation 
law (2.13) retains its form in the macroscopic description, i.e. that 



qt(f)-Uf) = fduj u {Vf) Wt.se R, / G V m (Q). 

J s 



(2.20) 
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3. Connection between the Quantum Picture, the Phenomenological Dynamics 

and the Hydrodynamical Fluctuations 

We now seek an inter-relationship between the quantum and hydrodynamical prop- 
erties of the macroscopic field qt(x) and its current jt(x) in the limit where iV tends to 
infinity. In order to formulate this limit, we shall henceforth indicate the iV-dependence of 
the quantum model by attaching the superscript (N) to the symbols E, 0, U, q,j, q 
and j. The symbol E, without that superscript, will be reserved for the limiting case where 
iV becomes infinite. The symbol Q, on the other hand, will continue to represent the fixed 
region occupied by E^, in the hydrodynamical scaling, for all N. 

Our basic assumptions concerning the relationship between the quantum and hydro- 
dynamic pictures of the model are that, in the limit iV^oo, 

(a) the stationary hydrodynamic fields q(x) and j(x) are the expectation values of the 
quantum fields q[ N \x) and j^ N \x), respectively, for the steady state 0^); and 

(b) the regressions of the fluctuations of these fields are governed, in a sense that will 
be made precise in Section 4, by the same dynamical laws (1.7) and (1.8) as the weak 
perturbations Sq t (x) and Sj t (x) of q(x) and j(x), respectively. 

The regression hypothesis (b) is a natural generalization of that proposed by Onsager 
[17] for fluctuations about equilibrium states. We remark here that, since V spaces are 
complete, these assumptions imply that the classical fields q(x), j(x), Sqt(x) and Sjt(x), 
introduced in Section 1, are distributions. 

3.1. Quantum Statistical Formulae the Hydrodynamical Variables. It follows 
immediately from our specifications that the above assumption (a) signifies that 

q(x) = hm^ 00 ($W,Qf ) (x)$W) (3.1) 

and 

j(x) = \™ N ^(*( N \jl N \x)* {N) ), (3-2) 

the t-independence of the r.h.s.'s of these formula being guaranteed by the stationarity of 
0^. 

In order to bring the hydrodynamical description of the model into line with thermo- 
dynamics, we introduce the field 9{x) = ($i(x), . ,9 m (x)), conjugate to q(x) as defined by 
the space-dependent version of Eq. (2.5), namely 

9(x) = s'(q(x)). (3.3) 

Since we are assuming that the system is perpetually in a single phase region, and thus 
that the function s' is invertible, it follows from this formula that the fields q(x) and 9(x) 
are in one-to-one correspondence. 

Turning now to the hydrodynamical equation (1.4), we see immediately that the sta- 
tionary field q(x) is determined by the requirement that J-{q; x) = 0, together with the 
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conditions imposed by the E^- 1 — 71 coupling at the boundary dfl of O. In order to specify 
these conditions, we denote by dVtj the section of dQ, where is in contact with TZj. 
We then assume the following boundary condition. 

(R) On the section dVtj of the boundary of E, the classical field 9(x) of this system 
takes the value 9j of the control variables of the equilibrium state in which TZj is initially 
prepared. Thus the array of reservoirs fixes the form of 9(x) and therefore of q(x) on <9fi. 

This assumption signifies that, on the hydrodynamic time scale and in the limit iV^oo, 
the local thermodynamical variables 9(x) of E spontaneously take up the same values as 
the reservoir with which this system is in contact at its boundary. The assumption is 
fulfilled by the models of Refs. [9-11]. 

Note on the Phenomenological Dynamics: V6> as Driving Force. In the general situa- 
tion where the field qt is time-dependent, we define its thermodynamical conjugate to be 
the field 9t given by the space-time dependent version of Eq. (2.5), namely 

9 t (x) = s'{q t (x)). (3.4) 

Thus, in view of our assumption that the system is perpetually in a single phase region, 
the function s' is invertible and the phenomenological law (1.4) may be expressed in the 
form 

^q t (x) = V.g(6 t ;x), (3.5) 

where the functional Q is determined by J according to the formula 

g{s'(q t );x) = -J(q t ;x). (3.6) 

In particular, in the case of nonlinear diffusion, it follows from Eqs. (1.4), (1.5), (2.5) and 
(2.6) that this phenomenological law reduces to the form 

-q t (x) + V.(K(9 t (x))V9 t (x)) =0, (3.7) 

where, in correspondence with the general relationship (2.5) between q and 6>, 

K{9) = K(q)J(q)^K([sr 1 (e)j([sr 1 (e). (3.8) 

One sees immediately from Eq. (3.7) that the gradient of the thermodynamical field 9 t 
acts as the hydrodynamical driving force. 

3.2. Linearized Perturbations of the Hydrodynamics. In view of our above 
remarks, 8q t is a distribution that satisfies Eq. (1.7) and vanishes on <9fi. We assume 
that the linear operator C appearing in that equation is the generator of a one-parameter 
semigroup, {T t \teR + } := T(R+), of transformations of £>' m (fi). The solution of Eq. (1.7) 
is then 

5q t = T t _Ms V t>s>0. (3.9) 
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Correspondingly, by Eq. (1.8), 

Sj t = K8q t = K,T t _Ms V t>s. (3.10) 

Further, by Eq. (3.9) and the dissipativity condition stated in the paragraph before Eq. 
(1-7), 

V - limt^ooTtV = V i)eV' m (Q) (3.11) 

or equivalently 

V-lim t ^ oo T t *f = 0V feV m ({l), (3.12) 

where {T t *|teR + } is the one-parameter semigroup of transformations of £> m (fi) dual to 
T(R + ). We denote its generator by £*, which is just the dual of C 

3.3. The Hydrodynamical Fluctuation Fields. We define the quantum fields, 
Zt N) ( x ) = (CiJ ( x ),~ >Cm,t( x )) and Vt = {ViP( x )>~ ^mj( x )), representing the fluctu- 
ations of the hydrodynamically scaled field q^ N \x) and the associated current j[ N \x), by 
the formulae 

e\ x ) = LT[^\ x )-{^ N) A N \ x )^ N) )l (3-13) 

and 

4 N) (x) = Lf [j <"> ( x ) - (*<"> , f t N) (x)^)} , (3.14) 

the normalization factor L^ 2 being natural for this scaling. The corresponding smeared 

fields £,f N \f) an d Vt N \9) are then the observables obtained by integrating these fields 
against test functions / (E.V m {Q)) and g (eVy(Q)), respectively. Thus, it follows from 

Eqs. (2.20), (3.13) and (3.14) that satisfies the local conservation law 

£t (JV) Cf) ~6 N Hf) = f ^^(V/jVMG R, / e V m (Q N ). (3.15) 

J s 



The dynamical properties of the fluctuation field Q are encoded in the correlation 
functions 

^(/ (1) ,,/ W ;ti,,t r ) = (* iN) ,^\f (1) ). 4?\f (r) )* (N) )- (3-16) 

This formula, together with Eqs. (2.15) and (3.13), serves to express in terms of the 

smeared fields q~t N \f) of Section 2. Thus, in view of our stipulation there that the actions 
on $( Ar ) of the monomials in these fields are continuous in the /'s, and £'s, it follows that 
W( N ' is continuous in all its arguments. Further, it follows from the stationarity of the 
state <f>( N ) and the self-adjointness of the observables ^ N \f) that 

.,/ W ;ti + a,. .,t r + a) = W^(f^\. .J^;t u . .,f P ) V aeR, (3.17) 



and 



W™(fW,. ,/< r >;fi,. ,t r ) = ^(/W. ,fi); (3.18) 
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while the positivity of <j)( N ^ implies that (A$>( N \ A^ N ^)>0 for any polynomial A in the 
smeared fields £ { t N) (f). Thus choosing A = ELi c fcdfi(/ (M) )- -St^ (f (k ' rk) ), where the 
c's are complex constants and p is finite, 

YZ, l=1 w wiN) V ik,rh) >- >/ (M) / (M) ,- ,/ (, ' r,) ;tfc,r fc ,- ,V,)>0- (3-19) 

3.4. Hydrodynamic Limit of the Fluctuation Process. We now assume that 
W( N > converges to a functional W in the hydrodynamic limit where iV^oo, i.e. that 

hm^o^CfW. .,f P ) = .,/ (r) ;ti,. ,t r ) 

V/ (1) ,.,/ (r) 6r(0), t l5 .,t r GR, rGN. (3.20) 

Hence, in view of the continuity properties of and the completeness of V spaces, W 

is continuous in the /'s and measurable in the t's. It is therefore a zero order distribution 
with respect to the latter variables [33]. Further, it follows immediately from Eq. (3.20) 
that W inherits the stationarity, Hermiticity and positivity properties of W^ N \ as given 
by Eqs. (3.17)-(3.19). Thus 

WtfM, . ., / (r) ; h + a, . ., t r + a) = W(f^\ . ., /«; t u . ., * r ) V aeR, (3.21) 

TF(/«. .,tr) = W(fW,. (3.22) 

and 

J*, ,ftCi» r (/ ( ' f ' r|) ) . ./^J^ 5 ,. l / (,,r,) ;^,. .tM^i-i, - »^)>0- (3-23) 
'fc,t=i 

It follows from these properties that, by Wightman's reconstruction theorem [25], 
corresponds precisely to the quadruple (H, V, £, where 

(a) H is a Hilbert space, 

(b) V is a unitary representation of R in H such that V t , the image of t (gR) under V, is 
strongly measurable; 

(c) £t(x) is a Hermitian operator valued distribution, of class V' m (Q), in H, which imple- 
ments the time translations of £, i.e. 

£ t+s (:r) = VtU^Vf 1 ; (3.24) 

and 

(d) \& is a vector in H that is invariant under Vt and cyclic with respect to the polynomials 
in the smeared fields £t(f) obtained by integrating £t(x) against £> m (n)-class test functions 
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The functional W is then related to the smeared field an d the cyclic vector \1/ by the 
formula 

W{p\. ,t r ) = ^^ tl (f {1) )- -6.(/ (r) )^)- (3.25) 

3.5. Conditions for W to represent a Classical Stochastic Process. The ques- 
tion of whether W represents a classical stochastic process reduces to those of whether 
(a) it defines a quantum stochastic process in the sense of Ref. [34] and (b) this pro- 
cess has the abelian properties of a classical one. Now the condition (a) is fulfilled if 
the smeared Hermitian fields &(/) are self-adjoint since, in this case, the unitary op- 
erators {exp(i£t(/)) \ fET> m (Q)} generate a VF*-algebra Aft and the correlation functions 
{(v&j.Ftj. .Ft r ^)\Ft s EAft s ', s = 1, . , r} define a quantum stochastic process, as formulated 
in [34] . Further, the classicality condition- 7 (b) is simply that of the intercommutativity of 
the operators £t(/). 

The following proposition provides a sufficient condition for the functional W to rep- 
resent a quantum stochastic process. 

Proposition 3.1. The functional W uniquely defines a quantum stochastic process 
indexed by D m (0)xR, if there is a bounded, positive functional (f,t)—>F t (f) on that 
product space such that 

\W(fU,. ,/< r >;fi,. ,t r )\<r 2 F tl (fV). .F u (f (r) ) V f (1 \ . ., /< r > e2> ra (fl); ^ . .,t r GR. 

(3.26) 

Comment. We shall subsequently establish in Prop. 6.1 that, under the assumptions 
of our scheme, the process £ is Gaussian. Since that result implies that the truncated r- 
point functions induced by W all vanish and thus that Eq. (3.26) is satisfied, it signifies a 
consistency of our assumptions. 

Proof of Prop. 3.1. As noted above, W defines a stochastic process if the Hermitian 
operators £t(/) are self-adjoint; and by Nelson's theorem [35], a sufficient condition for this 
is that each of these fields has a dense domain of analytic vectors. To prove that this is the 
case, subject to the assumption of Eq. (3.26), we note that it follows from that inequality 
and Eq. (3.25) that, for arbitrary /, /W, . , in £> m (fi) and t, t u . ,t r in R, 

\Umu(f (1) )- •6,.(/ (r) )^ll<(P + 2 F t (/)^ 1 (/ (1) )...^(/ M ) 

and therefore that the H-valued function z(eC)^J2^L z P M) p ^t 1 -6,.(/ (r) )^/P ! 
has an infinite radius of convergence. Hence, in view of the cyclicity of \P with respect to 
the polynomials in the smeared fields {&(/)}, these fields are self-adjoint and therefore W 
corresponds to a stochastic process. 

We shall assume henceforth that W does indeed define a stochastic process. In order 
to formulate a condition for its classicality, we introduce the following definition. 



J Here we consider classical processes as special (abelian) cases of the quantum ones. 
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Definition 3.2. (1) We define V (resp. to be the set of polynomials in the 

smeared fields {&(/) (resp. £ t (JV) (/)) |/e£> m (fi), teK} and we define the bijection P^pW 
of V onto by the prescription that p( N ) is the element of pW obtained by replacing 
£ by in the formula for P. 

(2) For PeP and iVeN, we define the vector } (eW^ w ) by the formula 

*W = pW*W (3 . 27) 

We now note that, by Eq. (3.25), the classicality condition that the operators £t(/) 
intercommute is equivalent to the invariance of W(f^\ . , f^;ti, . , t n ) under the permu- 
tations 

(/ w ,^)-(/ (r+1) ,Wi); 

and by Def. (3.2) and Eqs. (3.12), (3.16), (3.20), this latter condition may be expressed 
in the form 

lim^(*^,[^(/),^(/')]*g O )=0 

V PeV, /,/'eP m (fi), t,t'eR. 

Moreover, we can set t' = here without loss of generality, since is invariant under 
llj- N ^ and therefore, by Eq. (2.14), Def. 3.2 and the definition of ^ t N \ the manifold 
p( N )^( N ) is stable under this unitary transformation. Consequently, we have the following 
proposition, whose significance we shall discuss below. 

Proposition 3.3. Under the above assumptions, the process £ is classical if and only 
if £t(f) sa tisfi es the condition that 

lim^^, = VPeP, f,f'eV m (Q), teK. (3.28) 

Comment. In order to relate the condition (3.28) to the microscopic picture, we 
infer from Eqs. (2.10), (2.14)-(2.19) and (3.13) that this condition signifies the following. 

(1) In the case where t^O, 

hm^ooL^V™ / dxdy(^\ML N x,L%t),q l (L N y)}^ ) )f k (x)fl(y) = 

V f,f'eV m (n), PeV, (3.29) 

which is evidently a space-time asymptotic abelian condition on the field q. 

(2) In the case where t = 0, 

l™ N ^L- N ^^ p N \j( N \g fJ ,)^)=0 
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V f,f'ev m (n), Per, (3.30) 

where <?/,/' is the element of T>™(Q) whose rth component is 

9f,f';r = ^2 k hcrkifkVfl- (3-31) 

Thus, Eq. (3.30) signifies the avoidance of the catastrophe whereby, for fixed PeV, the 
expectation value of the smeared hydro dynamically scaled current p N '{gf,f) in the vector 
state would grow as rapidly as L 2 N with increasing N . 

4. The Stochastic Process £: Regression and Local Equilibrium Hypotheses 

and the Generalized Onsager Relations 

We now assume that the conditions of Props. 3.1 and 3.3 are fulfilled and hence that 
£ is a classical stochastic process, indexed by RxD m (fi). In a standard way, we denote 
the expectation functional of the random variables for this process by E. Thus, by Eq. 
(3.25), 

-M/ (r) )) = .&.(/ (p) )*) Vfi,. ,f P e R, e V m (Sl). 

(4.1) 

We note that, by Eqs. (3.20), (3.25) and (4.1), the process converges to £, i.e. its 
correlation functions converge to the corresponding ones for £, as N^oo. Further, in view 
of the observation following Eq. (3.20), the correlation function E^ tl (f^)- -^t r (f^ r ')) is 
continuous with respect to the /'s and measurable with respect to the t's. 

Conditional Expectations. For any random variable F of the ^-process and for teR, 
we denote the conditional expectations of F with respect to the cr-algebras generated by 
Mf)\feV m (n))} and {&(f)\t'<t, fEV m (n)} by E(F\£ t ) and E(F\£< t ), respectively. 

4.1. The Regression Hypothesis. This hypothesis is just the canonical gener- 
alization of that assumed by Onsager [17] for fluctuations about equilibrium states. Its 
essential import is that the evolution of a small hydro dynamical deviation from a steady 
state does not depend on whether the deviation has arisen from a spontaneous fluctuation 
or from a weak perturbation of the system fc . Thus, in mathematical terms, the regression 
hypothesis asserts that, for fixed s and t>s, the evolution of E(£t\£s) is governed by the 
same law as that of the linearised perturbation 5qt of the deterministic trajectory q t , i.e., 
by Eq. (3.9), that 

= [T t _ s £ s ](/)=£ s (Tt s /) V t>a. (4.2) 
Hence, since Nelson's forward time derivative [36] of &(/) is defined to be 

DUf) ■■= lim^ +0 «- 1 E(et +u (/) -WMt) (4.3) 

fc As in Onsager's theory, the assumption of this equivalence between the consequences 
of fluctuations and weak perturbations is not quite innocuous, since the modifications of 
the variables q due to the former are (3(iV _1 / 2 ), whereas those due to the latter are of 
order of a different small parameter that represents the strength of the perturbation. 



19 



and, since C is the generator of T(R+), it follows that 

Dit{f)=Cit{f). (4.4) 
Further, defining the static two-point function Ws ■ V m (Q) xD ra (0)^R by the formula 

W s (f, /') = £(£(/)£(/')) V /, fG^(O), (4.5) 
it follows from Eq. (4.2) and the stationarity of the £- process that 

E(tt(f )&(/')) = W s (T t *_ t ,f,f') V f,f'eV m (Q), t,t'(<t)eR. (4.6) 

4.2. Local Equilibrium Conditions. Our next assumption asserts essentially that, 
in a nonequilibrium steady state, the statistical properties of the fluctuation field £ in a 
'small' neighbourhood, Af(x), of an arbitrary point x (eO) simulate those enjoyed by these 
fields in the true equilibrium state corresponding to the value q(x) of the thermodynamic 
variable q. This is a mesoscopic local equilibrium condition, since it involves only the 
fluctuation field £ and is thus weaker than that of microscopic local equilibrium [37], 
which would signify that the microstate of £ in M{x) simulated the equilibrium microstate 
corresponding to q(x) there. Here we note that even this stronger condition has been shown 
to be fulfilled [38] by systems of fermions for which an Eulerian hydrodynamics has been 
established. Moreover, it may be expected to ensue more generally from the fact that 
the ratio of the hydrodynamic time scale to that of the microscopic processes (collisions 
etc.) is infinite, since that implies that local values of the hydro- thermodynamic variables 
q change negligibly in the time taken for the latter processes to generate equilibrium in 
macroscopically small spatial regions. 

In order to precisely specify our mesoscopic local equilibrium hypothesis, we start by 
formulating the relevant properties of hydrodynamical fluctuations about true equilibrium 
states for which the stationary classical field q(x) is assumed to be uniform. 

Equilibrium Fluctuations. We recall that, for a finite system, the equilibrium proba- 
bility distribution function, P, for macroscopic observables A is determined by the entropy 
S(A) according to the Einstein formula 

P(A) = const.exp(S(A)), 

and this serves to relate the static correlation functions for the fluctuations of these observ- 
ables to the thermodynamics of the system. The generalization of this relation to infinite 
systems has been derived by a quantum statistical treatment [15, Ch.7, Appendix C] of 
equilibrium states and takes the form 

where E eq is the equilibrium expectation functional for the fluctuation process, J(q) is 
defined by Eq. (2.6) and (., .) is the inner product on £> m (fi) defined by the formula 



(f,f') = ^ =1 J n dx f^fk( x )- 



(4.8) 
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It follows from Eqs. (4.2) and (4.7) that 

E eq {Uf)Uf')) = {T:_J,J(q)f) V fJ'E 2> m (fi), t> 8 . (4.9) 

Further, recalling the assumption, introduced in Section 2.5, of the invariance of the quan- 
tum field q( N \x) under the time-reversal antiautomorphism r and assuming that the 
equilibrium state z <f>iq of (E^) + 1Z) is likewise time-reversal invariant, it follows from 
the stationarity of this state and Eq. (3.13) that 

<^;*rV)e w (/')> = (*$ht™{f)&U)) = (^^ ( t N) (f')^ N Hf)). 

On passing to the limit of this equation as iV^oo, we see that 

and therefore, by Eq. (4.9), that 

E eq (Z(7tf)Z(f'))=E eq (Z(7tf')t(f), Vt>0. 
Consequently, since C* is the generator of T*(R + ), 

£ e<? (£(£*/)£(/')) =E eq ^(C*f)^f)) VfJ'eV m (Q). (4.10) 

Local Form of Equilibrium Correlations. We formulate the local properties of the 
equilibrium fluctuations in terms of test functions that are highly localised around an 
arbitrary point xq of Vt. Specifically, for f£.V m (VL), xq<E.VL and eGR+, we define the function 
f XOje on the Euclidean space X by the formula 

f X0 , e (x) = t- d ' 2 f{t-\x - x )) V x gO, feV m (n). (4.11) 

Since O is a bounded open subregion of X, it follows that the restriction of f XOj€ to O 
belongs to the space V m {VL) for sufficiently small e. In this case, we may take Eq. (4.11) to 
define a transformation /— >/ X0)€ of D m (fi), with e representing the degree of localization 
of the latter function about the point xq. 

We now note that, by Eqs (4.8) and (4.11), the r.h.s. of Eq. (4.7) is invariant under the 
transformation f—>f Xo ,e and therefore it follows from that equation that the equilibrium 
fluctuations enjoy the local property given by the formula 

iim ei0 ^(e(/ a; o,e)e(^ o ,e)) = (/, Aq)f) v x en, f, fev m (n). (4.12) 

Further, in the case of nonlinear diffusion, it follows from Eq. (1.10) that, for perturbations 
of the equilibrium state, C = K(q)A, with q constant. Hence, for fluctuations about 

1 The same assumption would not be valid for nonequilibrium states, since these gen- 
erally carry currents of odd parity with respect to time reversals. 
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equilibrium, it follows from Eq. (4.7) that both sides of Eq. (4.10) are invariant under the 
transformation f^f Xo , e , f'^f'x ,^ E eq —>e 2 E eq , and consequently 

\im el0 e 2 E eq (t(£*f Xo , e MfL o ,e)) = lim €io e 2 S efl (e(£*/i ,«)e(/xo, e ) ^ efl. (4.13) 

Local Equilibrium Conditions for Nonequilibrium Steady States. We now assume that, 
for these states, the natural counterparts of the local conditions (4.12) and (4.13) still hold, 
i.e. that 

\im el0 E(Z(f X0 MfLe)) = (/, A<lM)f) V x GO, /, f'eV m (Q) (4.14) 

and 

\im el0 e 2 E(a^f X0 , e )af X0 ,e)) =lim £l0 e 2 E(aC*f X0 ,M(f^)) V x GO, f,f'eV m (n). 

(4-15) 

These are our local equilibrium conditions, which manifestly concern the fluctuation field 
f only. 

4.3. Generalized Onsager Reciprocity Relations. The following proposition 
represents a generalization of the Onsager reciprocity relations to nonequilibrium steady 
states of the nonlinear diffusion process. 

Proposition 4.1. Under the assumption of the regression and local equilibrium hy- 
potheses, the transport coefficients of the nonlinear diffusion process satisfy the position- 
dependent Onsager relations 

K kl (6(x)) =Ki k (9(x)) VxgO, fc,/e[l,m]. (4.16) 

Proof. Since we employ the same argument as that for nonequilibrium states of 
conservative systems in Ref. [15, Ch. 7], we shall just sketch the proof here. We start by 
introducing the linear transformation £ of £> m (fi) by the formula 

C :=K(q(x ))A. (4.17) 

It then follows, after some manipulation, from Eqs. (1.10), (3.8), (4.14) and (4.17), together 
with the continuity properties of the functions K, J and q, that 

hm ei0 e 2 E(Z([£* - £$]/* ,e)£(O) = V /> f'eV m Q), x eQ. (4.18) 
This implies that C may be replaced by C in Eq. (4.15), i.e. that 

\im el0 e 2 E(^ f X0 , e )af X0 ,e)) =lim ei oe 2 i?(e(Ox ,J^(/xo, e )) V f,f'(=V m (Sl), x GO. 

(4.19) 

Further, since, by Eqs. (4.11) and (4.17), 

e 2 £(),/xo,e = [A)/]x ,e) 
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Eq. (4.19) reduces to the form 

It follows from this equation, together with Eqs. (3.8), (4.14) and (4.17) that 

(A/,K(£(* ))f ) = (A/',K(%))/), V f,f'eV m (Q), x en. (4.20) 
Further, since, by Eq. (4.8), 

(A/,/')=(A/',/) Vf,f<=V m (Q), 

and since the actions of A and if (#(:ro)) on £> m (fi) intercommute, Eq. (4.20) is equivalent 
to the following formula. 

(Af,K(9(x ))f) = (Af,K*(0(x o ))f) V/,/'GP m (0), x gO, (4.21) 

where K* is the adjoint of K. Hence, the matrix K{Q(xoj) is symmetric for all points xq 
in This is equivalent to the required result. 

5. Fluctuating Currents, Chaoticity and the Onsager- Machlup Process 

5.1. A Preliminary Observation. We now aim to extend the stochastic process 
£ so as to include the currents associated with these fluctuations. To this end we recall 
that, under the assumptions of Props. 3.1 and 3.3, £ t converges a classical process £, 
indexed by D' m (0)xR, with £t(/) continuous in / and measurable in t. We shall now 
argue that, by contrast, r]( N > cannot converge to a process r\ possessing the corresponding 
continuity and measurability properties. To show this, we suppose that the correlation 
functions for v^ N > converge to those of a process rj, indexed by Vy l (£l)x~R.. Then, since 
C* is the generator of T*(R+), it follows from Eqs. (3.15), (3.20), (4.1), (4.5) and (4.6) 
that t 

/ dsj ds 2 E( Vsi (Vf) Vs2 (Vf))=E(Uf)-af)] 2 ) = 
Jo Jo 

2£(£(/)[£(/)-£(T*/)]) = -2 f dsW s {f,Tt X* "/)) V JeV m (n), teK+. 

Jo 

Now the r.h.s. of this equation is 0(t), whereas the l.h.s. would be 0{t 2 ) if E (r} Sl (g)r} S2 (g)^ 
were continuous in g and measurable with respect to si and s 2 - Hence, we cannot assume 
that r]( N ) converges to a process r\ that possesses these continuity and measurability prop- 
erties. 

5.2. The Processes £ and rj. In view of this observation, we proceed somewhat 
differently, starting with the definition 

( ( t N s \g) ■■= fdu^\g) VgeVWl), t,seR. (5.1) 

J s 
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We assume that the cyclic vector <&( N ) lies in the domain of all monomials in the operators 
£u (/) and Q h a (g) as / and g run through V m {Vt) and T>y(Q), respectively, and t, s 
and u run through R. We further assume that the correlation functions given by the 
expectation values of these monomials for the vector state <f>( N > are continuous in their 
spatial test functions and time variables, that they converge pointwise to definite limits 
as iV-^oo, and that these limits satisfy the canonical counterparts to the assumptions of 
Props. (3.1) and (3.3). It then follows, by analogy with the arguments of Section 3, that 
the quantum process (^ N \ C^) converges to a classical one, (£, C), whose two components 
are indexed by D m (0)xR and I)™(0)xR 2 , respectively, and are continuous with respect 
to their spatial test functions and measurable with respect to their time variables. 

In view of Eq. (5.1) and the fact that the process ( is the limiting form of Q N > as 
iV^oo, we term C the time-integrated current. We note that since by Eq. (5.1), 

it follows that, correspondingly, 

Ct,s=Ct,u + Cu, s and Ct,t=0, (5.2) 
Further, by Eqs. (3.15) and (5.1), 

& (JV) (/) - 6 N) (f) = Cg°(V/) V /eP-(Q), f, sGR, 
and hence, correspondingly, 

W) - Uf) = Ct,.(v/) v fev™(n), t, sen, (5.3) 

which is just the local conservation law for £. 

5.3. Extension of the Regression Hypothesis: Secular and Stochastic Cur- 
rents. By Eq. (1.8), the increment djt in the phenomenological current due to a per- 
turbation Sq t of the field q t is )C5q t . Correspondingly, by way of extending the regression 
hypothesis of Section 3, we designate the secular part of the time-integrated fluctuation 
current (t,s to be 

C t s , e s c := f*dvK4 w (5.4) 

J s 

where /C, defined formally by Eq. (1.8), may now be interpreted as a mapping from V' m {Vt) 
into Vy 1 ^). We define the time- integrated stochastic current to be the residual part of 
Ct, a , namely 

Cm = Cm — Ct,s i 

i.e., by Eq. (5.4), 

Cm = Cm ~ / dulC£ u . (5.5) 

J s 
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In view of this formula, we may re-express the local conservation law (5.3) in the form 

W)-Uf) = fduU^Vf) + Cm(V/), 

J s 

or equivalent ly, since Eqs. (1.9) and (3.15) imply that VX = — £, 

W)~Uf)= f duZ u (C*f) + w t , a (f) V/GD m (0), t,seR, (5.6) 

J s 

where 

w t ,s(f) ■■= Ct,.(V/)) V /eX> ra (n), f, sGR. (5.7) 
Further, since, by Eqs. (5.2) and (5.7), 

w t ,s=w t ,u + w U)S and w t , t =0, (5.8) 

Eq. (5.6) is formally a Langevin equation. However, the condition for it to qualify as a 
bona fide Langevin equation is that w has the temporal stochastic properties of a Wiener 
process. The following proposition, which we shall prove in Appendix A, establishes that 
its two point function does have the requisite properties. Further assumptions concerning 
the chaoticity of the time-integrated stochastic current which will be introduced in 
Section 5.4, then lead to a picture in which w is indeed a fully fledged Wiener process. 

Proposition 5.1. Assuming the regression hypothesis, the local conservation law 
(5.3) and the definition ofwt, 

E(w tiS (f)Uf)) = V t>s>u, /, f GT?™(Q) (5.9) 

and 

E(w t M)^t',s'(f'))=-[Ws(^fJ')+W s (fX*f'Ms,t]n[s\t']\ 

V t, s(<t), t\ s'(<t')eR, /, f GP m (Q), (5.10) 

where the last factor represents the length of the intersection of the intervals [s,t] and 
[s',t'] and Ws is the two- point function defined by Eq. (4-5). Further the process w is 
non- trivial, i.e. w tjS does not vanish. 

5.4. The Chaoticity and Temporal Continuity Hypothesss. We assume that 
the stochastic current is chaotic in the sense that the space-time correlations of Ct,s(x) 
are of short range on the microscopic scale. This assumption is designed to represent 
Boltzmann's hypothesis of molecular chaos, as transferred from the local particle velocities 
to the stochastic currents. Since tends to infinity with N, it signifies that the space-time 
correlations of Ct,s(x) are of zero range on the hydrodynamic scale. Further, in accordance 
with the central limit theorem for fluctuation fields with short range spatial correlations 
[39] , we assume that the process ( is Gaussian. Thus, our chaoticity hypothesis is that 

(C.l) The process ( is Gaussian; 



25 



(C.2) E(C t , a (g)Ct>,A9')) = 0X(s,t)n(s',l/) = fr, and 
(C.3) EiCtMCt'Ag 1 )) = if supp((7)nsupp((7') = 0. 

It follows immediately from (C.l) that the process ( is completely determined by its 
two-point function E (( t , s (g)(t' ,s' (g')) ■ In view of the discussion following Eq. (5.1), this is 
continuous with respect to the test functions g and g' and measurable with respect to the 
time variables £, s, t' and s'. We now strengthen this conclusion by the following continuity 
hypothesis to the effect that it is continuous with respect to the time variables. 

(C) The two-point function E (( t ,s(g)(t' , s ' (g')) is continuous with respect to the time vari- 
ables £, s, s'. 

The following proposition, which we shall prove in Appendix B, stems from a appli- 
cation of a key theorem of Schwartz [33, Theorem 35] to the process £, subject to the 
assumptions (C.2) and (C). 

Proposition 5.2. Under the assumption of the hypotheses (C.2), (C.3) and (C), 
together with the condition of continuity with respect to its spatial test functions, the two- 
point function for the process ( takes the form 

E&Mtt'Ag')) =r(g,g')\[s,t}n[s',t'}\ Vg,g'eVy(n), t, s,t' , s'eR, (5.11) 

where reV'™(n)®V'™(Q) and suppress, x')eQ 2 \x' = x}. 

5.5. A Local Equilibrium Condition for the Currents. 

In order to extend our local equilibrium condition to the stochastic currents of the 
nonlinear diffusion process, we start by formulating the two point function at equilibrium 
for the process (. 

Equilibrium Two Point Function for (. Assuming again that the field q is uniform at 
equilibrium, we infer from Eqs. (1.6) and (1.7) that in this situation C = K(q)A, with 
q constant. Hence, by Eqs. (3.8), (4.7), (5.7) and (5.10), together with the symmetry of 
J(g), which follows from Eq. (2.6), 

E eq {(tAVf)(t>,s>(Vf')) =-[{Af,K(e)f) + {K(e)f,Af)]\[s,tns',t% 

which, by Eq. (4.16), is equivalent to the following formula for the unsmeared two-point 
function for (. 

d 2 

7 E eq (Ct, S ;kA x )Ct',s';lA^)) = -2K kl {e)AS{x-x')\[s,t}r)[s',t% (5.12) 



dx^dx u 



where Ct,s-k,^ is the //'th spatial component of the /c'th component of the field ( t , s = 
(Ct,s;i) •• 7 (t,s;m) and the summation convention is employed for the indices \x and v. Recall- 
ing now our assumption, at the start of Section 2, that the interactions are translationally 
and rotationally invariant, we assume that the corresponding symmetries are unbroken 
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in the pure equilibrium phase and thus that the process £ is invariant under the space 
translations and rotations that are implemented within the confines of O. We remark here 
that the limitation in Euclidean symmetry imposed by the boundedness of O is not serious 
from the physical standpoint, since O is an open subset of X and so any point of it, as 
viewed in the microscopic picture, is infinitely far from the boundary of S. 

Assuming then that the equilibrium two-point function for £ is invariant under space 
translations and rotations, we may express it in the form 

E eq ((t, S] kA x )(t',s'-,iAx')) = S k i{x-x')8^\[s,t)r)[s',t% (5.13) 

where SkiEV'(Q). It follows from this formula that Eq. (5.12) reduces to the following 
differential equation for Ski- 

AS kl (x) = 2K kl (9)AS(x). (5.14) 

Further, by condition (C.3) and Eq. (5.13), the distribution Ski has support at the origin, 
and therefore [33, Theorem 35] Ski(x — x', t) is a finite linear combination of S(x — x') and 
its derivatives. Hence the only admissible solution of Eq. (5.14) is 

S k i(x) = 2Kki(e)8(x) 

and therefore, by Eq. (5.13), the equilibrium two-point function for £ is given by the 
formula 

x — x')5^\[s,t]n[s' ,t']\ . (5.15) 
Equivalently, the equilibrium two-point function for the smeared field (t,s(d) takes the form 

E eq faMCf,.' (</)) = 2 {g, K{e)g') v \ [s, t]n[s\ t'} \ 

Vg,g'eV™(n), t,s,t',s'eR, (5.16) 
where (.)y is the inner product in T>y(Q) defined by the formula 

(d,9')v = YZ-i I dxg(x).g'(x)V g,g'evy(n)= 

^ — y m ^— yd f 

Z^ fc=1 Z^ M=1 J n dx 9kA x )9kA x )^ ( 5 - 17 ) 

and where gk,^ is the /u'th spatial component of gk- 

Local Property of the Equilibrium Two Point Function. We formulate the local prop- 
erties of the stochastic current £ along the lines employed in Section 4.2 for the process £. 
Thus, for (xo, e)GOxR + , and e sufficiently small, we define the transformation g^g X() ,e of 
Vy(Cl) by the formula 

9x M = e- d ' 2 g{e-\x - x )). (5.18) 
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We then observe that, by Eqs. (5.17) and (5.18), the transformations t^e 2 t, g^g Xo ,e, of 
the times and test functions lead to the multiplication of the smeared two-point function 
of Eq. (5.16) by the factor e 2 . Thus, 

e- 2 E eq (C eH ^ s (g Xo , e )C(g' Xo , e )) = 2{g,K(6)g') v \[s,t]n[s\t']\ 

Vx gO, g,g'eV^(Q), t,s,t',s'eK. (5.19) 

The local property of the two-point function for ( at the point xq is then obtained by 
passing to the limiting form of this equation as e^O. 

Local Equilibrium Property for the Stochastic Current in the Nonequilibrium Steady 
State. In view of the last observation, we assume that, in the nonequilibrium steady state, 
the process £ enjoys the local equilibrium property obtained by passing to the limit e^O 
and replacing E eq and 6 by E and 9(xq), respectively, in Eq. (5.19). Thus we assume that 

\im e ^ e- 2 E(C eH , e 2 8 (9x ,e)Cen',e^(g' Xo ,e)) = 2(g , K (e(x ))g') v \[s, t]n[s' , t']\ 

Vx gO, g,g'eV^(n), t, s(<t), t', s'(<t')eR. (5.20) 

This is our local equilibrium condition for the stochastic current. 

5.6. Explicit Form of the Two Point Function for (. By Prop. 5.2, this function 
is determined by the functional T, which by Eqs. (5.11) and (5.20), possesses the following 
local equilibrium property. 

hm el0 T(g X0 ^g' X(he ) =2(g,K(e(x ))g') v y x en, g,g'eV^(n). (5.21) 

The following proposition, which will be proved in Appendix C, provides an explicit formula 
for the functional T, which stems from a combination of the chaoticity condition (C.3) 
and the local equilibrium condition (5.21). 

Proposition 5.3. Under the previous assumptions, together with the local equilibrium 
condition of (5.21), Y is given by the formula 

T(g, g') = 2(g, K e g') v V g, g'eV^(Q), (5.22) 

where Kg is the matrix-valued operator KoQ in Vy{Vt), i.e. 

K e (x)=K(6(x)). (5.23) 

The following corollary is an immediate consequence of this proposition and Prop. 

5.2. 

Corollary 5.4. Under the same assumptions, the two-point function of the stationary 
process ( is given by the formula 

E(CUg)Ct',Ag'))=2(g,Kog')v\[s 1 t]n[s\t']\ 
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V g, g'eV^(Q), t, s(<t), t', S '(<OeR- (5.24) 



5.7. The Generalized Onsager-Machlup Process £. 

It now follows immediately from Cor. 5.4 and Eq. (5.7) that 

E(w t , a (f)v>t>Af')) = 2(Vf,KoVf') v \[s,t]n[s',t'}\ 

V /, f'eV m (Q), t, s, t', s'eR. (5.25) 

Hence, by the chaotic hypothesis (C.l) and Eqs. (5.7) and (5.25), w is a generalized 
Wiener process. Further, on re-expressing Eq. (5.6) in the form 

dZt = £Ztdt + dw t ,s, (5.26) 

we see that, in view of the additive property (5.8) of w, the fluctuation field £ executes 
a it generalized Onsager- Machlup process; while Eq. (5.25) signifies that the two-point 
function for w corresponds precisely to that assumed for the stochastic force in Landau's 
fluctuation hydrodynamics [18]. 

In order to derive the properties of the process £ from those of w, we note that, since 
C is the generator of T(R + ), the solution of the Langevin equation (5.26) is given by the 
formula 

& = T t _ s £ s + I T t _ u dw u , s V t, s(<t)eR, (5.27) 

J s 

or equivalent ly, 

&(/) = ZsW-sf) + [ dw U)S {Tt_J) V /GI?™(Q), T, s(<t)eK+. (5.28) 

The following proposition, which we shall prove in Appendix D, is a natural general- 
ization of standard properties of the Brownian motion of a single particle that ensue from 
the Langevin equation governing its velocity (cf. [36]). 

Proposition 5.5. Under the above assumptions, 

(1) £ is a Gaussian Markov process, and 

(2) the fields Wt,s and ^ u are statistically independent of one another if s and t are greater 
than or equal to u. 

Comment. It follows from this proposition and Eqs. (4.5) and (4.6) that the process 
£ is completely determined by the forms of the semigroup T*(R + ) and the distribution 

w s . 



6. Long Range Spatial Correlations of the £- Process 
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6.1. The Static Two-Point Function for £. By Eq. (4.5), the unsmeared form of 
the D /m (Q)<g>'D /m (fi)-class distribution Ws is given by the formula 

W s (x,x') = E(Z(x)®Z(x')). (6.1) 



The following Proposition provides an explicit formula for Ws, as well as a differential 
equation for this distribution in terms of the semigroup T*(R+), and the transport function 
K e . 

Proposition 6.1. Under the above assumptions, 

/•OO 

W s (f, /') = 2 / dt(VT*f, K e VTtf) v V/, f'eV m ({l) (6.2) 
Jo 

and, further, the generalized function W~s(x,x') satisfies the equation 

[£®I + I®C']W s (x, x') = 2V. (Ke(x)VS(x - x')) , (6.3) 
where C is the version of C that acts on functions of x' . 

Proof. By Eq. (4.5) and the stationarity of the ^-process, 

w s (f, /') = E(UfMf')) v /, f'ev™(n), teR + 

and therefore, by Eq. (5.26), 

Ws(fJ') = E(t(T*f)t(T*f')) + [* E(t(T*f)dw Vi0 (T*_ u f')) + 

Jo 

1 E(Z(T;f)dw u ,o(T:_J))+ f f E(dw u , (T t *_ u f)dw u ,, (T:_ u ,f')) 

Jo Jo 

V/,/'Gl? m (0), tEK+. (6.4) 

Now, by the dissipativity condition (3.12), the first term on the r.h.s. of this equation 
vanishes in the limit t^oo, while by Eq. (5.9), the second and third terms there vanish. 
Hence, it follows from Eq. (6.4) that 

W s (f, /') = linw, /' f E{dw ufi {Tt_ u f)dw u ,, {T;_ u ,f')) V /, f'eV m (Q). (6.5) 
Jo Jo 

Further, by Eq. (5.25), 

E(dw u , (f)dw u ,,o(f')) = 2{Vf,K e Vf') v 8{u - u')dudu' 
and consequently Eq. (6.5) reduces to the form 

W s (f, f') = limt^oo2 / du(VT t *_J,VT t *_J') v =2 I du(VT*f, V2^/')v, 

Jo Jo 



L 
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which is equivalent to the required formula (6.2). 

Further, since £* is the generator of T*(R+), it follows from Eq. (6.2) that 

Ws(£*f,f') + Ws(f,£*f') = 2 J dt-(VT t *f,K e VT t *f') v 

and consequently, by the dissipativity condition (3.12), 

W s (C*f, /') + Wstf, C*f) = -2(V/, K e Vf')v V /, f'eV m (n), 
which, by Eq. (6.1), is equivalent to the required formula (6.3). 

6.2. Long Range Spatial Correlations. In order to provide a precise character- 
ization of long range correlations, we first recall that the ratio of the macroscopic length 
scale to the microscopic one is infinite. Consequently, correlations of finite range on the 
microscopic scale are of zero range on the macroscopic one. Accordingly, we term the range 
of correlations 'short' or 'long' according to whether or not it reduces to zero in the macro- 
scopic picture. Thus our condition for long range spatial correlations for the £-field is simply 
that the support of the distribution Ws does not lie in the domain {(x, x')eQ 2 \x = x'}. 
The following proposition establishes that the spatial correlations of £ for the nonlinear 
diffusion process are generically of long range. 

Proposition 6.2. Let <& q be the m-by-m matrix-valued function on Q defined by the 
formula 

^ q (x) = AK e (x) + V.^ q (x), (6.6) 

where 

Then under the above assumptions, a sufficient condition for the spatial correlations of £ 
to be of long range is that either $ 9 does not vanish or that the matrix ^ q is symmetric. 

Comments. (1) The Proposition establishes that the correlations are generically of 
long range, since the specified conditions on $ 9 and ^> q can be satified only for special 
relationships between the functions Koq and soq; and these are generally independent of 
one another, since s and K govern the equilibrium and transport properties, respectively, 
of E. By contrast, the corresponding correlations for equilibrium states are generically of 
short range, except at critical points. A treatment of critical equilibrium correlations of 
fluctuation observables is provided by Ref. [40]. 

(2) In the particular case of the symmetric exclusion process [9-11], n = 1, d = 
1, K(q) = 1, s(q) = —qlnq — (1 — q)\n(l — q) and q(x) = a + b.x, where a and b (^0) are 
constants. Thus, in this case, it follows from Eqs. (1.6), (2.6), (6.6) and (6.7) that ^ q = 
and $q(x) = — 26 2 ^0. Hence, long range correlations prevail in this model, in accordance 
with the results obtained by its explicit solution in Refs. [9-11]. 



o 

.dan (x) 



Kkk'(q{x)) [Jin(q(x))Vq k '{x) - J k 'i(q{x))Vqi>(x)]. (6.7) 
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Proof of Prop. 6.2. Suppose that the static spatial correlations of £ are not of long 
range, i.e. that the support of the distribution Ws lies in the domain {(x, x')eQ 2 \x' = 
x}. Then it follows from this supposition and the local equilibrium condition (4.14), by 
precise analogy of the derivation of Eq. (5.24) from corresponding conditions of zero range 
correlations and local equilibrium for the process £, that 

W s (x, x) = J q (x)5(x - x), (6.8) 

where 

J q {x) :=J(q(x)). (6.9) 
Hence, by Eqs. (1.10), (3.8) and (6.7)-(6.9), 

(C®I)W(x, x') = A[K e (x)5(x - x')} + V.[-V q (x)5(x - x')]. (6.10) 

Further, by Eq. (6.1), 

(I®£')W(x,x') = [(£'®I)W(x',x)] tr , 

where the superscript tr denotes transpose, and therefore, by Eq. (6.10), 

(I®C')W(x, x) = A'[K 9 (x' - x)8(x' - x)] tr + V'.[%(x')6(x' - x)] tr ', (6.11) 

where A' and V are the versions of A and V, respectively, that act on functions of x' . 
Consequently, since Kg is symmetric, by Eqs. (4.16) and (5.23), it follows from Eqs. (6.6), 
(6.10) and (6.11) that 

[C®I + I®C']W s (x,x') = 
2V.(K (x)VS(x - x')) + $ q (x)5(x - x') + [V q (x) - ^ q r (x)].VS(x - x'). (6.12) 
On comparing this equation with Eq. (6.3), we see that 

$ q (x)5(x - x') + [$ q (x) - ^l r (x)].V5(x - x') = 0, 

i.e. that $ ? vanishes and that ^ q is symmetric. These, then, are conditions that ensue 
from the assumption of short range correlations of the ^-process. We conclude, therefore, 
that the violation of either of these conditions signifies that the correlations are of long 
range. 

7. Concluding Remarks. 

We have proposed a macrostatistical treatment of nonequilibrium stady states of quan- 
tum systems that is centred on the fluctuations of their hydrodynamical variables. The 
key physical assumptions on which this treatment is based are 

(a) the regression hypothesis for the hydrodynamic fluctuation field £; 

(b) the chaoticity of the associated currents, as represented by their time integrals (t,s', 
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(c) the local equilibrium conditions on the stochastic process comprising £ and £; 

(d) the space-time scale invariance of the phenomenological equation of motion (1.4), as 
exemplified by the case of nonlinear diffusions; and 

(e) the invariance of the quantum field q, and correspondingly of the classical field £, under 
time reversals. 

On the basis of these assumptions and certain technical ones, we have obtained a picture 
that provides natural generalizations of the Onsager reciprocity relations and the Onsager- 
Machlup fluctuation process to nonequilibrium steady states, together with a demonstra- 
tion that the spatial correlations of the hydro dynamical variables are generically of long 
range in these states. Furthermore this picture is expressed exclusively in terms of the 
phenomenological functions representing the equilibrium entropy, s(q), the transport coef- 
ficients K{6) and the hydro dynamical boundary conditions. This may easily be seen from 
the comment at the end of Section 5, together with Eqs. (1.10), (3.8) and (6.2) and the 
fact that the semigroup T(R) is completely determined by its generator C. 

Let us now discuss the assumptions (a)-(e) a little further. In our view, for reasons 
expressed in Sections 4.1, 4.2 and 5.4, the first three of these seem natural from the physical 
standpoint, though they are very hard to prove in concrete cases. On the other hand, it 
is clear that assumptions (d) and (e) are not universally valid: for example, they both fail 
in the important case of Navier-Stokes hydrodynamics. Consequently, it is of interest to 
consider how the macrostatistical picture presented here might be extended to situations 
where (d) and (e) are replaced by weaker assumptions. 

In fact, the weakening of (e) provides no serious problems, since the locally conserved 
fields of continuum mechanics are generally either even or odd with respect to time reversals 
[41]. Accordingly, we replace (e) by the assumption that each of the quantum fields qj has 
either even or odd parity with respect to time reversals, i.e. that 

Tqj(x) = Rjqj(x), Rj = ±1, j = l,. .,ra, (7.1) 

where again r is the time-reversal antiautomorphism. This weakened assumption then 
leads to the nonlinear version of Casimir's extension [41] of Onsager's theory, wherein Eq. 
(4.16) is modified to the formula 

K kl (9(x )) = R k RiK lk (R6(x )), (7.2) 

where 

R(6) := .,R n 6 n ). (7.3) 

Similarly, the modification of assumption (e) to the form given by Eq. (7.1) presents no 
serious problems for the other issues treated here. 

On the other hand, there does not appear to be any natural generalisation of the scal- 
ing assumption (d), which lay behind the interdependence of the ratios of the macroscopic 
to microscopic scales for distance and time, the former ratio being and the latter L 2 N 
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(or more generally Lj^). Moreover, one sees from Eqs. (2.15) and (3.13) that this inter- 
dependence was essential to the limit procedures of Eqs. (3.1) and (3.20). Nevertheless it 
does not appear to be essential to the key physical ideas that 

(i) the ratios of the macroscopic to microscopic scales for both distance and time are 
extremely large, and 

(ii) the currents associated with the locally conserved quantum fields satisfy the chaoticity 
assumption of Section 5.4, whereby the space-time correlations of their fluctuations decay 
within microscopic distances and times. 

Since such chaoticity does not necessarily require any interdependence of the ratios of the 
macroscopic to microscopic scales for distance and time, it appears reasonable to expect 
that some version of the present macrostatistical model should still be applicable even in 
the absence of macroscopic space-time scale invariance. 

Thus, from the standpoint of mathematical physics, a most challenging question is 
whether the present scheme can be generalized to a setting which does not require the 
scale invariance of the macroscopic law (1.4). Presumably such a generalization would 
require a difficult multi-scale analysis. 



Appendix A: Proof of Proposition 5.1 

We shall first prove Eqs. (5.9) and (5.10) and then demonstrate the nontriviality of 
the process w. 

Since £ is the generator of T(R + ), Eq. (5.9) follows immediately from Eqs. (4.2) and 
(5.6). 

It then follows from Eqs. (5.6) and (5.9) that the l.h.s. of Eq. (5.10) vanishes if the 
intervals [s,t] and [s',t f ] do not intersect. Hence, in view of Eq. (5.8), the proof of Eq. 
(5.10) reduces to that of the same formula with s = s' and t = t' and t>s. Thus it suffices 
for us to prove that 

E(w t M)wtAf')) = - {Ws(C*f, /') + W s (f, £*/')) I* " *l 

VM(<*)eR, fj'ev m (n). (A.l) 

We start by inferring from Eq. (5.6) that the l.h.s. of Eq. (A.l) is the sum of the 
following four terms:- 

E[(Uf)-Uf))(Uf')-Uf'))}, (a) 

- f duE[{Uf) ~ Uf))U^ f% (6) 

J s 

- f duE[UC* f){W') - Uf))] (c) 

J s 
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and 

ft pt 



j du j dvE(U£*f)U£*f))- (d) 

J s J s 

Since t>s and the ^-process is stationary, it follows from Eqs. (4.5) and (4.6) that 

Term (a) = 2W s (f, /') - W s (T?_ a f, /') - W s (f, T?_J% (A.2) 

Term (b) = - f duW s (T t *_J, £*/') + f duW s (f, T*_ s £*/'), (A3) 

J s J s 

Term (c) = - f duW s (£* f,T t *_J') + f ' duW s {T* u _ a C* f,f) (A A) 

J s J s 



and 



r-t rU pt ft 

Term(rf)=/ du dvW s (T* u _ v C j \ C f) + du dvW s (C*f,T*_ u C*f). (A5) 

J s J s J s J u 

Since W5 is linear in each of its arguments and since C* is the generator of T*(R + ), 
it follows that Eqs. (A. 3-5) may be re-expressed in the following forms. 

Term (b) = - f duW s (T t *_J, C*f) + W s (f, T t *_ s f) - W s (f, /'), (A.6) 

J s 

Term (c) = - f duW s (C*f, T t *_J') + W s (T t *_J, /') - W s (f, f) (A.7) 

J s 

and 

Term (d) = f du [-W s (f, C*f) + W S {T* U _J, C* f) + W s (C*f, T t *_J') - W 3 (£*f, /')] • 

J s 

(A.8) 

It follows now from Eqs. (A.2) and (A. 6-8) that the sum of the terms (a), (6), (c) and 
(<i), which comprises the l.h.s. of Eq. (A.l), is equal to the r.h.s. of that equation. This 
completes the proof of Eq. (A.l) and thus of Eq. (5.10). 

Finally, we employ a reductio ad absurdum method to establish the nontriviality of 
the process w. Thus, we assume that w tjS vanishes. It then follows from Eq. (5.7) that 

W s (C*f, /') + W s (f, C*f) = V /, f'eV m (Q) 

and hence that 

W s (C*T?f,T?f') + Ws(T?f,C*T?f') = V /, fel? m (Q), teR+. 
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Since Ws is linear in each of its arguments and since £* is the generator of T*(R+), this 
signifies that 

jW s (T t *f,T t *f')=0 

and therefore, since T = /, that 

w s {TtfXf) = Ws(fj') v f,fev m (n), ten+. (A9) 

Moreover, by Eq. (4.5) and the dissipativity condition (3.12), the l.h.s. of Eq. (A. 9) 
vanishes in the limit t^oo. Hence Eq. (A. 9) implies that the static two-point function 
Ws vanishes. This conflicts with the fact that, by Eqs. (4.5), (4.11) and (4.14), 

\im 6l0 Ws(f w f^ £ ) = {f,J(q(x ))f), 

which does not vanish identically. This contradiction establishes that the assumption of 
the triviality of w is untenable and thus completes the proof of the proposition. 

Appendix B: Proof of Proposition 5.2. 

We start by noting that, in view of Eq. (5.2) and condition (C.2), the proof of this 
proposition reduces to that of the formula (5.11) for the particular case where s = s', t = t' 
and s<t. Thus we need only prove that 

E(( t , s (g)( t , a (g'))=T(g,g')(t-s)V g,g'eV™(n), t,s(<t)eR, (B.l) 

where V is an element of V™®V™ with support in the domain {{x, i')6fi 2 |x' = x}. 
To this end, we start by defining 

F g , g ,(t,s):=E(( t MttA9')) (B.2) 
and inferring from Eq. (5.2) and condition (C.2) that 

F g , g ,(t, s) = F 9 , g ,(t, u) + F g79 ,(u, s) for t>u>s. (B.S) 

Further, by Eq. (B.2) and the stationarity of the process C, 

F gj9 ,(t, s) = F gj9 ,(t + b,s + b)V beK, 

which signifies that F g ^ g i may be expressed in the form 

F g , g ,(t, S ) = F g , g ,(t-s) Vs,teR, (BA) 

where, by condition (C), F g ^> is a continuous function on R. It follows now from Eqs. 
(B.3) and (B.4) that 

F 9 , 9 '(t) + F 9,At') = F 9,9'(t + V t, t'eR + (B.5) 
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and hence that 

F g ,g'(nt) = nF g , g/ (t), V teR+, neN 

or equivalently 

F g ,g,(t) =riFg t g,{t/ri), VtGR+, fl'eN\{0}. 

These last two equations imply that 

F 9 ,Art)=rF g , g ,(t) 

for all non-negative t and positive rational r; and further, by condition (C), this result 
extends to all positive r. Hence the action of F g ^ g > on R + takes the form 

F g , g ,(t) =T(g,g')tVteR+, (B.6) 

where T(g,g') := F g , g ,(l). By Eqs. (B.2) and (B.4), Eq. (B.6) is equivalent to the 
required formula (B.l); and, moreover, it follows from condition (C.2) and the continuity 
and linearity of the l.h.s. of that equation with respect to the test functions g and g' that 
T is indeed an element of V'™®V™ with support in the domain {(x, x')eQ 2 \x' = x}. 

Appendix C: Proof of Proposition 5.3. 

We base the proof of Prop. 5.2 on the following lemma. 

Lemma C.l Let Q\ be any open subset offl whose boundary, dQi, does not intersect 
dVt. Then, under the assumptions of Prop. 5.2, the restriction of the two-point function 
T to the spatial domain Vt\ is given by a finite sum of the following form. 

V 9l! ,'el>!?(«i), t,s,t',s'eR, (C.l) 

where 

(i) the C's are continuous functions on O with support in some arbitrary neighbourhood of 

n i; 

(H) 9k, (j, is the fx'th spatial component of the k'th component of g = (g\, . ,g m ); and 
(Hi) for n= (m,. .,n d )eN d , := d ni+ -+ nd ldx\\. dx n d d . 

Proof of Prop. 5.3 assuming Lemma C.l. We start by inferring from Eq. (5.18) 
that, for any g,g'&T>y(Q), xo&Q and e sufficiently small, one can find an open subset Q\ 
of O such that g X(jy6 and g' X(j lie in Vy (Hi). Hence, by Eqs. (5.18) and (C.l), 

r(<7x ,e> 9x ,e) = 
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Vg,g'eVV(n). (C.2) 

where \n + n'\ := Ylt=i( nk + n 'k) : evidently the effective domain of integration here is 
supp((7)flsupp((7 / ). Since the functions C are continuous, the summand on the r.h.s. of 

this equation will diverge, as e^O, unless either n and n' are both zero or C k \ u (xo) = 0. 
Hence the local equilibrium condition (5.21) implies that the only non-vanishing C's are 
those for which n and n' are zero. Thus, Eq. (C.l) reduces to the form 

rWHEw-X" -i / ^SM^^^-M^K.^v^^'e^ini). (C.3) 

Correspondingly, Eq. (C.2) reduces to the form 

r(^ ,e, g' X0 , e ) = Y^ l=1 J2l „ =1 / dxC kLA Xo + ex )9kAxHA x ) 

Vg,g'eT>V(a). (C.4) 

It now follows immediately from this formula and the local equilibrium condition (5.21) 
that 

Em ^ — / / / / 

,L / dxC k]l^A x °)9kA x )9iA x ) = 2{g,K(9(x ))g ) v 

Further, in view of Eq. (5.17), this last equation signifies that 

<t^>) = 2 M*(*))<w (c.5) 

and consequently that Eq. (C.3) reduces to the required formula (5.22), at least for 
g, g'eVy^x). The extension to all g, g' in Vy (O) is trivial, since for any pair of elements 
of the latter space, one can always choose fii to be an open subset of that space that 
contains their supports. 

Proof of Lemma C.l. Since the test functions gk,^ and g[ v in Eq. (C.l) are 
arbitrary elements of T>(Q) : this lemma reduces to the following one. 

Lemma C.2. Let T be a V (fi 2 )- class distribution whose support lies in the region 
{(x, x')e.Vt 2 \x' = x} and let fii be an open subset of fi whose boundary, dfli, does not 
intersect <9fi. Then the restriction of T to the domain {f®f'\f,f'^T>{Q,i)} is given by a 
finite sum of the form 

r{f®f) = E „ Nd / dxC^'(x)d:f(x)d:'f'(x) V fJ'eV^), (C.6) 

^— ' n,n' E!>i Jq 



where the C's are continuous functions on fi with supports in some neighbourhood of Q.\. 
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Proof of Lemma C.2. Let a be a 22(0)- class function which takes the value unity 
in fii and whose support lies in a compact connected subset, K, of O whose boundary, 
dK, does not intersect either 90 or 90i. We define the distribution T (gD^O 2 )) by the 
formula 

f(x, x') = a(x)a(x')T(x, x'). (C.7) 
Thus, T coincides with T in Of and 

su PP (t)c{(x, aZ^V = x}. (C.8) 

We define $ to be the linear transformation of X 2 given by the formula 

$(y,z) = (y + z,y-z)Vy,zeX (C.9) 

from which its follows that 

$ _1 (x, x') = (-(x + x'), ^(x - x')) V x, x'el. (CIO) 

We then define 

6 := ^(O 2 ) = {(y,z)EX 2 \(y±z)en}, 
and we define the bijection F^F of £>(0 2 ) onto V(@) by the formula F = Fo$, i.e. 

F(y, z) = F(y + z,y-z)V (y, z)eQ. (C7.ll) 

Correspondingly we define the distribution T (gX>'(G)) in terms of T by the formula 

f(F) = f(F) V Fg£>(0 2 ). (C.12) 

It follows from Eqs. (C.8), (C.ll) and (C.12) that 

supp(t)cKx{0}. (C.13) 

We want to restrict T to an open subset of G which contains the support of this 
distribution and takes the form 2 x J, where 2 and J are open subsets of O and X re- 
spectively. Accordingly, we choose b to be a positive number that is less than dist(dK, 90), 
the minimal distance between the boundaries, dK and 90, of K and O. We then define 
O2 := {y£X\(y, z)gG V |z|<&} and J := {zGX||z| < b}. It follows from these definitions 
that O2 and 2 x J are open subsets of O and G, respectively, that _ftTc0 2 and that 90 2 , 
the boundary of 2 , does not intersect either dK or 90. Hence, by Eq. (C.13), 2 x J is 
an open neighbourhood of supp(T) and the restriction, T', of T to this domain carries all 
the information we require. It follows from its definition that T / Gr ,/ (0 2 x J). 

Now let e be an arbitrary element of £>(0 2 ). Then for e'eV(J), T' induces a contin- 
uous linear functional 7^ on V( J) according to the formula 



fl(e') = f'(e®e') Ve'eP(J), 



(C.14). 
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where the mapping e— of V(Q,2) into T>'(J) is continuous. Further, it follows from Eqs. 
(C.13) and (C.14) that TJ has support at the origin and consequently, by Schwartz's point 
support theorem [33, Theorem 35], that this distribution is a finite sum of derivatives of 
5(z), with coefficients given by linear continuous functionals of e, i.e. 

tV) = E T»(e)(0V)(O), (C.15) 

where each T n eV(J). Further, in view of the definition of %!, it follows from Eqs. (C.13) 
and (C.14) that T n has support in the compact K and therefore, by Schwartz's compact 
support theorem [33, Theorem 26], it is a finite sum of derivatives of continuous functions 
on f^2 with support in an arbitrary neighbourhood of K. Consequently, by Eq. (C.15), 
the action of T' on X>(Ox J) is given by a finite sum of the form 

f '( P ) = J2, I dyD n '' n (y)d^d:F(y,z) z=0 \/ FeV(n 2 xJ), (C.16) 

where the D's are continuous functions on fi 2 with support in a neighbourhood of K. 
Hence, as T' is just the restriction of T to T>(Q,2X J) and since T coincides with T in Of, 
it follows from Eqs. (C.9)-(C.12) that Eq. (C.16) is equivalent to the formula 

T ( F ) = E dxC^'(x)d^;F( Xl x% /=x y fJ'eV(^), (C.17) 

which in turn is equivalent to the required Eq. (C.6). 



Appendix D: Proof of Proposition 5.5 
Part (a). The characteristic functional for the process £ is 

CCfW. .,tr) = S[exp(t^ =1 U(/ ( * ) ))] 

V/W,..,/MeD m (n); ,t r GR, reN. (D.l) 
Equivalently, since the process is stationary, 

C(/« . .,/<••>;*!,. .,*,) = S[exp(i^ =i et fc+ to(/ (fc) ))] 

V/^. .,/^Gl? m (0); t ,*i,. -,^eR, rGN. (£>.2) 

Here we are at liberty to choose to to be any real number and, for any specified set of 
times ti, . ., t r , we choose it so that t± + 1 , . , t r + t are all positive. It then follows from 
Eq. (5.28) that 

rtk+to 

(f (k) ) = m k+ t f (k) ) + / dw u , (Tl +t0 _jW) 

Jo 
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and therefore that Eq. (D.2) may be re-expressed as 

C(/W ..,/W;ti,..,t P ) = 

We now define 

C(fW,. .,/ (r) ;to,ti,- -,tr) =exp(i^ fc=i jf ^u,o(T* +to _ u /W)) (D.4) 

and, using the Schwartz inequality, we infer from the last two equations that 
|C(/«. ,/W;ti,. .,tr)-C{f<- 1 \..,f< r hto,t 1 ,. .,t r )| 2 

<^[|exp(^l =1 ^^ + *o/ (fc) ))-l| 2 ] 

^[(£l =1 emw (fc) )) 2 ] = E; il=1 £ («w w )f(W ))- 

It follows from the dissipativity condition (3.12) that the r.h.s. of this estimate vanishes 
in the limit to^oo, and therefore that 

C(/W,. ti,. = liin t0 ^ 0O C(/< 1 >,. .,/W;t ,ti,. .,*r), 

i.e., by Eq. (D.4), that 

C(/«. .,f P ) = lnn t0 ^ OOJ B[ex P (i^ fc=i y ^,o(T/ fc+to _ u /^))] . 

Since, by Eq. (5.7) and the chaoticity condition (C.l), the process iy is Gaussian, it follows 
immediately from this last equation that the process £ is Gaussian. 

In order to show that it is also Markovian, we need just to prove that, for teR and any 
random variable B> t generated by {£t t (/)|/££ ,m (^), u>t} : the conditional expectations 
of B> t with respect to the random variables for time t and for times <t are equal, i.e. that 

E(B>t\Z t ) = E(B> t \S< t ). (D.5) 

Now the random variables over the times >, = and <t are generated by linear combinations 
of terms 

F> t = exp(^ =i e ufe (/( fe ))), (D.6) 
F t =exp(z6(/)) (D.7) 

and 

F<t = exp(i^ =1 6, (/ /(0 )), (^-8) 



41 



respectively, where Uk>t>si and f^ k \ f and /'^ are elements of P m (0). 

It follows from Eqs. (D.l) and (D.6)-(D.8), together with the Gaussian property of £, 

that 

£^>t^) = • (D.9) 

and that 

E(F> t F< t ) = C(f l \ ..,/<">;«!,. .,« p )C7(/'< 1 ',. Sl , . .,s r )x 

Further, since Uk>t>si, it follows from Eqs. (4.5) and (4.6) that the summands appear- 
ing in the exponents in Eqs. (D.9) and (D.10) are equal to (2£ fc _ t / (fc) )f (/)) and 
i?(£(T* fc _ S; /( fc ))£(/ /(7 ))), respectively, and therefore those equations may be re-expressed 
as 

E(F> t F t ) = C(/W,. .,/<">;«!,. .,^)C(/;t)exp[-^ =i £;(e(r4_ t /W)e(/))] (D.ll) 
and 

E(F>,F< t ) = C(/«. .,/<">;«!,. .,u p )C(/'W,. .,/';<>„. .,s r )x 

-p[-El =1 EL 1 B (« T «-- fW )«- f ' ( ' > »]- 

Further, since E(F>t\£t) is the unique random variable of the ^-process at time t for which 

E(E{F> t %)F t ) =E(F> t F t ) 

for all F> t and F t of the forms given by Eqs. (D.6) and (D.7), respectively, it follows from 
Eq. (D.9), together with the stationarity and the Gaussian property of the process, that 

S(Fa,l6) = CTC,_,/a»...,r,;.,/. P »;o....o) eitp (E^^ (r ">-' / ' (fll3) 

Hence, by Eq. (D.8), 

E(E(F> t \£ t )F< t ) = 

Further, in view of the Gaussian property of the process, the last factor in this formula is 
equal to 

C(/'W, . ., /'<'•> » r )C(El =1 r i-</ W ; ')>< 

-p[-£lXL 1 B te( T »V t / ( ' , x»,(/' ( ' ) ))]- 
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Therefore since, by Eq. (4.6) and the semigroup property of T*(R,+), the summand in the 
exponent in this expression is equal to E^(T* k _ Si f^ k ^(f'^)), it follows from Eqs. (D.IO) 
and (D.14) that 

E(E(F> t \£ t )F< t ) = E(F> t F< t ). 

Hence 

E(F> t \S t ) = E(F> t \Z< t ), 
which signifies that the process is temporally Markovian. 

Part (b). Since Eq. (5.8) implies that wt, s = —Ws,t and since wt yS and £ u are 
Gaussian random fields whose means are zero, it follows from Eq. (5.9) that the latter two 
fields are statistically independent of one another if s and t are both greater than or equal 
to u. 
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